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Abstract 

We report our results on non-periodic experimental time series of pressure in a 
single cylinder spark ignition engine. The experiments were performed for different 
levels of loading. We estimate the noise level in internal pressure calculating the 
coarse-grained entropy from variations of maximal pressures in successive cycles. 
The results show that the dynamics of the combustion is a nonlinear multidimen- 
sional process mediated by noise. Our results show that so defined level of noise in 
internal pressure is not monotonous function of loading. 



Key words: Spark ignition engine, combustion, noise estimation 



1 Fax: +48-815250808; E-mail: g.litak@pollub.pl 
Preprint submitted to Elsevier Science 



8 February 2008 



Torque control system 



Cylinder pressure sensor 




Pressure acquisition 



Fig. 1. Schematics of the engine pressure measurement and control setup. 



1 Introduction 



The cycle-to-cycle combustion variability has been a subject of interest for 
many years [1,2,3,4,5,6,7,8,9,10,11,12]. It is possible that many different dis- 
turbances influence the process making it stochastic but in case of high sensi- 
tivity on the process conditions one should also consider appearance of chaotic 
behaviour [2,3,4,7,8,9,10,11]. The pressure signal can be used to control the 
engine combustion [15]. Its variations might originate from a complex dy- 
namics leading, presumably with some stochastic component, to non-periodic 
behaviour. Thus, the crucial problem is to understand the nonlinear dynam- 
ics of the process observing the internal pressure inside the cylinder. In the 
present note we discuss experimental results of a direct measurement of in- 
ternal pressure and analyse its time series using the coarse grained entropy 
method [13]. The method has been already used successively by Kaminski et 
al. [12] to examine the dynamics of combustion in a four cylinder combustion 
engine to examine the influence of an spark advance angle on cycle-to-cycle 
variations of heat release. In this paper we present analysis of pressure con- 
cerned to a single cylinder engine. Experimental works have been performed 
in the Engine Laboratory at University of Trieste. 

The paper is organised as follows. After the present introduction we will discuss 
experimental facilities in Sec. 2. Following Sec. 3 is devoted to analysis of 
pressure time series. In this section we show the experimental results and 
perform their initial analysis. In Sec. 4 we present the methodology of noise 
level estimation and we finally apply the coarse-grained entropy method to 
internal pressure time series. At the end we provide conclusions in the last 
Sec. 5. 



2 Experimental Facilities 



Tests were carried out on a SI - four stroke, single cylinder, Aprilia/Rotax 
engine. The instrumentation which was used during our test is schematically 
presented in Fig. 1. 
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The pressure sensor Kistler 6053C, was preferred for its small dimensions. 
Mounting the pressure sensor was a complex task because of the cylinder head 
profile. The sensor was placed in the relatively small available space, which 
remains between the five valves areas. It was centrally positioned, having the 
active part very closed to the spark plug electrodes. As the upper end of 
the sensor, stays under the cylinder head cover, a protective tube for the 
connecting cable was mounted. 

The encoder (Lika C58L), due to the engine configuration, was fixed on the 
camshaft. 

This solution generates complications regarding signal processing. Because 
of the transmission chain, there is no rigid connection between the engine 
crankshaft and the encoder shaft, therefore synchronisation problems may 
appear. To overcome this, before each testing session, motored engine tests 
were accomplished, to determine the shifting between the encoder trigger and 
pressure maximum peak. 

The reference level for the cylinder pressure was established using the intake 
pressure signal, acquired with an inductive pressure sensor. Before each test 
session, this sensor was calibrated using a simple system containing: a piston 
pump which creates the calibration pressure, a high precision manometer, a 
signal amplifier which amplifies the sensor output and a voltmeter [14]. 



3 Time Series Analysis 

Using our experimental standing we have measured internal pressure time se- 
ries for 1000 cycles measuring 6000 points per cycle. In Fig. 2 we have plotted 
pressure against the crank angle <fi for first five consecutive cycles of measured 
engine work against the crank angle for different torque loadings F = 0, 20 
and 40 Nm. Note that with larger loads the peak pressure is higher and the 
combustion zone concentrates near the TDC. Variations of pressure reflects 
the efficiency of combustion process in successive cycles. Such a phenomenon 
is associated with the burned fuel fluctuations and influence negatively on its 
consumption. The cyclic fluctuations of pressure in combustion process can 
be described by its maximal value P max in each cycle. Interestingly, as it was 
proposed in Ref. [15], it is possible to use the maximal values of pressure 
pmax i nvo i vec i j n control process to achieve stable combustion. The method 
of peak-to-peak analysis appeared to be useful to analyse some kind of dy- 
namical systems [16], as it is in our case. This peak-to-peak variation is our 
starting point for further analysis by means of a Coarse-Grained Entropy and 
an autocorrelation function. 
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Fig. 2. Pressure in first five cycles of measured engine work against the crank angle 
<fi for different torque loadings F = 0, 20 and 40 Nm at 4000 rpm. Numbers in 
figures 1-5 correspond to pressure curves in successive cycles, respectively. 
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Table 1 

Definitions of variables and symbols used in the paper. 



number of measured pressure consecutive in averaging procedure 


r 


internal cylinder pressure 


Pi 


maximum pressure in cycle 


pmax p 
r i ui 1 max 


average of maximum pressure series 


Pmax 


maximum pressure vector in the embedding space 


-pmax 

i 


Heaviside step function 


e(z) 


embedding time delay in cycles 


m 


cycle number 


hi 


embedding dimension 


n 


number of considered points in time series 


TV 


loading torque 


F 


crank angle 


4> G [0, 720°] 


threshold 


£ 


correlation integral 




coarse-grained correlation integral 


C n {e) 


correlation entropy 


K 2 


coarse-grained correlation entropy 


K2(e) 


calculated from time series coarse-grained entropy 


K 

J1 noisy 


correlation dimension 


D 2 


Noise-to-Signal ratio 


NTS 


standard deviation of data 


&DATA 


error function 


Erf(z) 


fitting parameters 


k , a, b 


standard deviation of noise 


a 


top dead center 


TDC 
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Fig. 3. Cycle-to-cycle variations of maximal pressure values P ma x of 900 cycles at 
different loads F = 0, 10, 20, 28, 40, 43 Nm, respectively. Note different scales in 
lower and upper panels. 




Fig. 4. Standard deviation a data of Pmax signals versus an applied torque F. Note, 
the arrow at F = 20 Nm indicates the first sudden grow of a data for relatively low 
loading. 
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To show its cycle-to-cycle changes we have plotted P max against successive 
cycles % in Fig. 3. One should note that both, average value of pressure P ma x 
as well as its standard deviation (Jdata (Fig. 4) increase with growing torque 
F. Interestingly for about F = 20 we observe the sudden growth of &data- 
This is a result of some diversity of a maximal pressure origin as pressure itself 
is a complex product of compression and combustion phenomena. Clearly 
one can see the intermittent behaviour of P ma x (Fig. 3 for F = 20 and 28 
Nm) can be explained in these terms. Lean combustion where the maximum 
of pressure is produced largely by piston compression with possible ignition 
missing corresponds to lower level of P max (Fig. 3) while effective combustion 
happens from time to time corresponds to spikes in Fig. 3. Similar effects can 
be seen for F = 28 and F = 40 but not for F = 43 where the distribution of 
Pmax is going back to a Gaussian type. 



4 Noise Level Estimation 



In the n dimensional embedding space [17] the state is determined by a mul- 
tidimensional vector 

r>max r pmax pmax pmax pmax \ /i\ 

1 — \ r i ! r i+m! f i+2m)-) r !+(n-l)mJ) \ L ) 



where m denotes the embedding delay in terms of cycles. The correlation 
integral calculated in the embedding space can be defined as [18,19] 

i N N 

c n ( £ ) = ^EE^- \\ p T ax - Pfl l)> (2) 



where N is the number of considered points corresponding to pressure peaks 
in cycles and is the Heaviside step function. For simplicity we use maximum 
norm. The correlation integral C n (e) is related to the correlation entropy K2(e) 
and the system correlation dimension D 2 by the following formula [18,19] 

lim C n (e) = D 2 \ne-nmK 2 (e). (3) 



The coarse-grained correlation entropy can be now be calculated as 
V( \ v 1 Cn ( £ ) dlnC"(5) 

K 2 (e) = hm In - , ;/. sa — — . (4) 

v 1 «^oo C n+1 (e) dn v 1 



In such a case the correlation entropy is defined in the limit of a small threshold 
e. 
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Fig. 5. Level of noise estimated for different loads (Torque). 



In presence of noise described by the standard deviation a of P™ ax time series, 
the observed coarse-grained entropy K noisy [13] can be written as 



K n0 i S y (e) 



m 



x 1 + 



£ 

2^ 



J hie + [k + 61n(l - ae)] 
e 2 /3 + 2a 2 -e/Vs' 



(5) 



Function g(z), present in the above formula, reads 

2 ze~ z2 
g(«) = - 



7T Erf (z) 



(6) 



where Erf(.) is the Error Function. The parameters k, a, b as well as o are 
unconstrained. They should be fitted in Eq. (5) to mimic the observed noisy 
entropy calculated from the experimental data. 

Using the above procedure in analysing the time series of pressure peaks we 
have calculated the coarse-grained correlation entropy for different loading 
levels (F = 10, 20, 28, 40, 43Nm). Fitting it to the last formula (Eq. 5) it was 
possible to estimate the noise level a. In Fig. 5 we present the obtained results 
with the help of NTS parameter defined as follows 



NTS 



a 



&DATA 



(7) 
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6. Autocorrelation function of the maximal pressure signal. 



500 



where odata is the standard deviation of maximum pressure signal data (Fig. 
4). The error on Fig. 5 comes from the fitting procedure of Eq.(5) to ob- 
served coarse-grained correlation entropy. Note, the leading tendency is that 
for higher loading the noise is larger, what is expected for most of the engines 
with lean combustion. However the cases of F = 20 and 43Nm could be the 
exception of this rule. At F = 20 we simultaneously observe the first sudden 
increase of <r data- It is natural that the overall effect (Eq. 7) will cause de- 
creasing of NTS. Interestingly, the nodal or very low loading F = signal (of 
Pmax) appears to be influenced more by a compression phenomenon mediated 
by combustion process (Fig. 2a) than combustion itself. In this case P max is 
highly correlated. In Fig. 6 we show the autocorrelation function for different 
torque loading F = 0, 20 and 40Nm. Note that for the lowest loading (F = 0) 
the autocorrelation function does not decrease in contrast to the other cases. 



On the other hand at F = 43 Nm (Fig. 5) we have to do with a different 
situation. In that case due to increasing loading the fresh fuel feeding is of 
higher level and cycles which misses of ignition do not happen any more. On 
the other hand for F = 40 and 43 Nm we observe further growing of the signal 
square deviation (Tdata with an increasing load. In such conditions it is not a 
surprise that in such a case we also observe a lower level of NTS. 
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5 Summary and Conclusions 

We applied the Coarse- Grained Entropy Method for the experimental data 
estimating the noise-level. This method enable to distinguish determined sig- 
nal, including chaos with its short time prediction scale, and random noise. 
Analysing our experimental data we decided to use maximal values of pres- 
sure P max . In all considered cases we have found a relatively high noise level. 
However it has appeared that larger loading led to more noisy system. For 
larger power the engine needs more fuel. Increase of fluctuation in combustion 
can happen if the engine work is optimised for the lowest fuel consumption 
or the exhaust gas toxicity. Another interesting point, discussed in the paper, 
was connected with a missing ignition effect. In fact such effects are present 
for small loading (Figs. 2a-b) but influence strongly P ma x distribution only if 
loading is large enough. 

We have also found that the pressure signal P max was more correlated (Fig. 
6) for small loading. Note, the noise level obtained for a correlated signal case 
by our method can be underestimated by a few percents. A noise component 
for small loading can be also related to a well known problem of idle speed 
engine control difficulties [20,21]. 
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